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Abstract 



Particle production in the forward region of heavy-ion collisions is shown to be due 
to parton recombination without shower partons. The regeneration of soft partons due 
Q^ ' to momentum degradation through the nuclear medium is considered. The degree of 

degradation is determined by fitting the p/p ratio. The data at ^fs = 62.4 GeV and 
.^j- . rj = 3.2 from BRAHMS on the px distribution of average charged particles are well 

\ reproduced. Large proton-to-pion ratio is predicted. The particles produced at any px 

should have no associated particles above background to manifest any jet structure. 

o 

in ■ PACS numbers: 25.75.Dw 

o 
o 

1 Introduction 

In an earlier paper [I] we studied the problem of hadron production in the transfragmentation 
region (TFR) in heavy-ion collisions. It was stimulated by the data of PHOBOS [2] that 
show the detection of charged particles at rj' > 0, where rf — rj — ybeam- We broadly refer to 
the rj > region as TFR. However, since the transverse momenta pr of the particles were 
^ ■ not measured, it has not been possible to determine the corresponding values of Feynman 
x, in terms of which TFR can more precisely be defined as the region with x > 1. More 
recently, BRAHAMS has analyzed their forward production data at ^/s = 62.4 GeV with 
both rj and pt determined [3]. It is then possible to interpret BRAHMS data by applying 
the formalism developed in [TJ, which is done entirely in the framework of using momentum 
fractions instead of rj. In this paper we calculate the proton and pion distributions in x and 
Pt and conclude not only that the p/n ratio is large, but also that there should be no jet 
structure associated with the particles detected at any pt in the forward region. 

In pQ the x distributions of p and 7r have been calculated for 0.6 < x < 1.2 in the recom- 
bination model [HEJE], taking into account momentum degradation of particle constituents 
traversing nuclear matter [7] and the recombination of partons arising from different beam 
nucleons. However, we have not considered the regeneration of soft partons as a consequence 
of momentum degradation in the nuclear medium. Since such soft partons significantly in- 
crease the antiquark distribution in the mid- a; region, it is important to include them in the 
determination of the pion distribution. Furthermore, no consideration has been given in [1] 
to transverse momentum, which is the other major concern in this paper. 
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In the following we use forward production to refer specifically to hadrons produced at 
x > 0.3, with the fragmentation region (FR) being 0.3 < x < 1, and TFR being x > 1. Any 
hadron produced in the TFR cannot be due to the fragmentation of any parton because of 
momentum conservation, since no parton can have momentum fraction > 1, if we ignore 
the minor effect of Fermi motion of the nucleons in a nucleus. In the FR hadrons with any 
Px that are kinematically allowed can, in principle, arise from the fragmentation of hard 
partons; however, the momenta of those hard partons must be even higher than the detected 
hadrons in the FR, and the probability of hard scattering into the region near the kinematical 
boundary is severely suppressed [8] . Moreover, there is the additional suppression due to the 
fragmentation function from parton to hadron. Thus the fragmentation of partons at any px 
in the FR (despite the nomenclature that has its roots in reference to the fragmentation of 
the incident hadron) is highly unlikely, though not impossible. The issue to focus on is then 
to examine whether there can be any hadrons produced in the FR with any significant px- 
If so, then such hadrons at any px would not be due to fragmentation and would therefore 
not have any associated jet structure. 

In contrast to the double suppression discussed above in connection with fragmentation, 
recombination benefits from double support from two factors. One is the additivity of the 
parton momenta in hadronization, thus allowing the contributing partons to be at lower x 
where the density of partons is higher. The other is that those partons can arise from different 
forward-going nucleons, thus making possible the sum of their momentum fractions to vary 
smoothly across x = 1, thereby amalgamating FR with TFR. These are the two attributes 
of the recombination process that makes it particularly relevant for forward production. Its 
implementation, however, relies on two extensions of what has been considered in [TJ , namely, 
the regeneration of soft partons and the transverse-momentum aspect of the problem, before 
we can compare our results with BRAHMS data [3J. 

It is useful to outline the logical connections among the different parts of this work. 
First of all, the degree of degradation of forward momenta through the nuclear medium 
is unknown. The degradation parameter k can be determined phenomenologically if the x 
distributions of the forward proton and pion are known, but they are not. We have calculated 
the x distributions for k = 0.6 and 0.8 as typical values serving as benchmarks. Since the 
normalization of the px distribution, which is known from BRAHMS data [3J , depends on the 
x distribution, k can be determined by fitting the px distribution. However, what is known 
about the px distribution is only for all charged particles, not p or ir separately. If there 
were experimental data on the p/it ratio (which is not yet available for the 62.4 GeV data 
that have the px distribution), one could disentangle the species dependence. Fortunately, 
there exist preliminary data on p/p and K/tt ratios at 62.4 GeV. We shall therefore calculate 
the p and p distributions and adjust n to render the ratio Rp/ P to be in the vicinity of the 
observed ratio. We shall then show that the px distribution of all charged particles can be 
well reproduced in our calculation. 

We shall assume factorizability in px and px dependences. The two are treated in Sec. 
II and III, respectively. The main difference between Sec. II and the earlier work in pQ is 
the inclusion of the regeneration of soft partons, a subject we now address. 
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2 Regeneration of Soft Partons 



Let us first recall some basic equations from Ref. pQ, which we shall refer to as I. Equations 
I-(16) and I-(32) give the proton and pion distributions in x (with the pt variables integrated 
out) for AB collisions in the recombination model 



Hf(x) = ^^^(x^,^,,^^), (1) 

J Xi X 2 X3 

H* B (x) = f^^F q f(x 1 ,x 2 )R 7T (x 1 ,x 2 ,x). (2) 

where the hadronic x is 2pi / i/i and the partonic Xi are momentum fractions. The recombi- 
nation functions R p and R w are given in pQ. The partons are assumed to arise from different 
nucleons in the projectile nucleus A and thus contribute in factorizable form of F AB , i.e., 
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(x 1 ,x 2 ,x 3 ) = F^ Xl )F^x 2 )F^x 3 ) , (3) 
F q f(x 1 ,x 2 ) = Fg( Xl )Fg(x a ) , (4) 

where v is the average number of wounded nucleons that a nucleon encounters in traversing 
the nucleus B at a particular impact parameter, given in Eq. I-(49). The effect of momentum 
degradation on the parton distributions is contained in the expressions 
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FSfa) = 1^ dy'G' D (y')K , (5) 
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F*(xi) = dy'G' D (y')L' q ^ 1 , (6) 

00 -v 

G'Av') = E«" 2 ^(K-V)-Te- P , (7) 

where k is the average momentum fraction of a valon after each collision and G(y) is the 
valon distribution in momentum fraction y before collision [9j. K(z) and L' q (z) are the quark 
distributions in a valon, with 

K(z) = K NS (z)+L' q (z) , (8) 

Kns{z) being the valence-quark distribution and L'Jz) the saturated sea-quark distribution 
after gluon conversion. This briefly summarizes the essence of determining the x distributions 
of protons and pions produced in AB collisions. 

To describe how the above should be modified in order to take into account the regen- 
eration of soft partons, we need to fill in the steps on how sea-quark distribution L' q (z) is 
derived. In addition to the valence quark in a valon, there are also sea quarks (q), strange 
quark (s) and gluons (g), whose distributions are denoted by Li(z),i = q, s,g. Their second 
moments satisfy the sum rule for momentum conservation [H [9] 

K NS {2) + 2 \2L q {2) + L s (2)l + L g {2) = 1 . (9) 
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Gluon conversion to qq changes the sea-quark distribution to 

L ' q , s {z) = ZiL g ,s{ z ) > 



(10) 



whose second moment satisfies the modified version of Eq. ([9]) where L g (2) is absent, i.e., 



K NS (2) + 2[2L , q (2) + L' s (2) 
From these equations we can determine Z\, getting 

4(2) 
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Z x = \ + 



2L q (2) + L s (2) 



(11) 



(12) 



This is what we obtained and used in I to calculate the hadron distributions. 

The degradation effect is parametrized by k such that 1 — k is the fraction of momentum 
lost by a valon after a collision. After v collisions, the net momentum fraction lost is 1 — k v . 
That fraction is converted to soft partons so that the new sea-quark distributions L" a 
satisfy a sum rule that differs from Eq. 
for conversion, i.e., 



q,s\ 



TTT) by the addition of extra momentum available 



K NS (2) + 2 [2Z;'(2) + L" s {2)] = 1 + (1 - k v ) . 
Assuming that only the normalization is changed, we write 

which yields, upon using Eqs. (jHJ) and (|T3l) . 

-k" + L b {2) 



Z 2 (k, v) — 1 + 



2L q (2) + L s (2) 



(13) 



(14) 



(15) 



In [T] we have considered the cases: k = 0.6 and 0.8 for b = 1 fm (0-5%) and 8 fm (30-40%). 
For any given b, the average v is known [see Eq. I-(12), (13)]. The dependence of v on v is 
Poissonian, as expressed in the last factor in Eq. (|7j). We now replace L'(z) in Eqs. © and (jHJ) 
by L"(z, k, v) and obtain the new distributions Fp(xi, k) and Fp(xi, k) defined in Eqs. (E]) and 
(Q, in which the summation over v in G' p (y') is now extended to include the v dependence of 
L''(z, k, u). As an illustration of our results on the effects of degradation and regeneration, 
we show in Fig. 1 the it-quark, Fg(x), and w-antiquark, F£(x), distributions before and 
after regeneration for b = 1 mb and k = 0.6. Note that with or without regeneration all 
distributions are highly peaked at x = because momentum degradation pushes all valons 
to lower momenta by a factor of k v (which for v ~ 6 is ~ 1/20). Regeneration increases 
Fp (x) significantly for x < 0.3, as shown by the dashed-dotted line above the dotted line. 
For F£(x), because of the dominance of the valence quark distribution K^s{z), the increase 
is minimal, as the dashed line is nearly all covered by the solid line. Similar changes occur 
for the d and d distributions. 

In the same way as we have done in [1] we calculate the proton and pion distributions 
in x for k = 0.6 and 0.8 and for 6 = 1 and 8 fm. The results are shown in Figs. 2-5. Since 
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the regenerated soft partons do not affect the hadron distributions for x > 0.8 (remembering 
that the hadron x is the sum of the parton Xj), we have plotted these figures for the range 
0.3 < x < 0.9. We emphasize here that the large x behavior in the TFR is not the central 
issue in this paper any more, as it was in pQ. 

In Figs. 2-5, in addition to our present result with regeneration (solid and dashed lines) 
we show also our previous result obtained in [1] without regeneration for the case k = 0.6 for 
the purpose of seeing the effect of regeneration. Note that the proton distributions in Figs. 

2 and 3 are not affected very much by the regeneration effect, but the pion distributions in 
Figs. 4 and 5 are increased. At x — 0.6 the increase is roughly around a factor of 3. 

In [TJ there was no data to compare with the calculated result on the x distributions. In 
particular, the degree of momentum degradation was unknown. Now, BRAHMS data show 
the pt dependence at r\ = 3.2 [3]. In order to fit the px distributions of the hadrons produced, 
we must have the correct normalizations, which in turn depend on the x distributions that 
we have studied. 

3 Transverse Momentum Distribution 

Having determined the longitudinal part of the hadronic distributions above, modulo the 
value of k, we now proceed to the transverse part. We have treated the degradation and 
regeneration problems on rather general grounds without restricting the x values and with 
Pt integrated so that pt never appears in our consideration of the x distribution. It is then 
natural to make use of that result in a factorizable form for the inclusive distribution 



which is, of course, an assumed form that is sensible when there is negligible contribution 
from hard scattering. 

For the transverse part, Vh(pr), we follow the same type of consideration as developed in 
[6], where particle production at intermediate pt is shown to be dominated by the recombi- 
nation process. Similar work in that respect has also been done in (TUKU]. In the absence of 
hard scattering there are no shower partons. Without shower partons there are only thermal 
partons to recombine. Thus for pion production we have TT recombination, while for proton 
we have TTT recombination, where T represents the thermal parton distribution [6] 



In the above equation p iT is the transverse momentum of ith. parton; Ci and T are two 
parameters as yet undetermined for the forward region in Au+Au collisions. In view of the 
factorization in Eq. (fl6l) we use the term thermal in the sense of local thermal equilibrium of 
the partons in the co-moving frame of a fluid cell whose velocity in the cm system corresponds 
to the longitudinal momentum fraction x. However, the value of T can include radial flow 
effect. 




(16) 




(17) 
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Limiting ourselves to only the transverse component, the invariant distributions of pro- 
duced pion and proton due to thermal-parton recombination are 

dN th 

ocC q C 9 exp(-Pr/T), (18) 



dN th 

cxC 3 p T exp(-p T /T), (19) 



where the proportionality factors that depend on the recombination functions are given in 
[6]. At midrapidity, thermal and chemical equilibrium led us to assume C q = Cq, and we 
have been able to obtain p/n ratio in good agreement with the data. Now, in FR (and in 
TFR) we must abandon chemical equilibrium, since q cannot have the same density as q, 
when x is large. But we do retain thermal equilibrium within each species of partons to 
justify Eqs. (ITSl) and (fl9l for the pt dependence. We join the longitudinal and transverse 
parts of the problem by requiring 

C q cc F«( Xi , k), C q cc Fl( Xj , k), (20) 

where Fp(xi,K) and Fp(xj,n) are the quark and antiquark distributions in their respective 
momentum fractions already studied in Sec. II above. The proportionality factors in the 
two expressions above are the same. Equation (|20|) connects the parton density from the 
study of the longitudinal motion to the thermal distribution in the transverse motion. Sub- 
stituting those relations into Eqs. (TT81) and (TT9"j) . and letting Fjj(xi, k), F£(xj, k) and other 
multiplicative factors be absorbed in the formulas for Hh(x,n) developed in [I], we obtain 
for the transverse part of Eq. (fT6l) 

V 7T (p T )=clexp(-p T /T), (21) 

V p (p T )=c 3 p p T exp(-p T /T), (22) 

where and c p are two proportionality constants to be determined by the normalization 
condition 

/ dp T pt V h (p T ) = 1 , (23) 
Jo 



i.e., 



= l/T, c p = 1/(2 1 / 3 T). (24) 

It follows from Eqs. (Tl6|) and (|23|) that we recover the invariant ^-distribution 

dNh s , . 

x-r± = H h (x, «) 25 
dx 

without undetermined proportionality factors. 
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The exponential factors in Eqs. (12T1) and (122]) give the characteristic behavior of hadrons 
produced by the recombination of thermal partons [6j [8]. Such exponential behavior is 
overwhelmed by power-law behavior at intermediate pr due to thermal-shower recombination 
when x is small and when light quarks contribute to the hadrons produced. However, when 
x is large, the shower partons are absent due to the suppression of hard scattering, so the 
exponential pt dependence becomes the prevalent behavior in the forward region. Since the 
data of BRAHMS [3] exhibit the pt distribution for a narrow range of 7] around 3.2, we can 
readily check whether Eqs. ffl6|) . (12T1) and (1221) are in accord with the data. We consider 
only the most central collisions for which H n (x, k) and H p (x, k) have been calculated in the 
previous section for 6 = 1 fm. The data of the pt distribution in [3] are, however, given for 
average charged particle (h + + h~)/2. To be able to make comparison with that, we need 
information on the magnitudes of contributions from K and p. Preliminary data on the K/ir 
ratio is ~ 0.15 [12], and that on the p/p ratio is ~ 0.05 |13j . We shall use the former ratio 
and calculate the latter. 

The enhanced q distribution enables us to compute Hp(x, k) exactly as in Eq. (pQ), except 
that Fp(xj) in Eq. ([3]) is replaced by F£(xj). The p/p ratio is constant in p T , since both 
p and p have the same px dependence given in Eq. (T221) . The value of the ratio, however, 
depends on H p (x, k) and H p (x, n). The value of x for both p and p is chosen to be 0.55 for 
reasons to be explained below when we discuss the pr distribution. Since the data on p/p 
are preliminary and imprecise at this point, we consider two values of k and obtain 

k = 0.76, R p/p = 0.031, 

k = 0.72, i2p /p = 0.058. (26) 

These results bracket the observed value of p/p at ~ 0.05 for y/s = 62.4 GeV and r\ = 3.2 
rrS]. Note that with a 5% decrease in k there is over 80% increase in p/p. This is a direct 
consequence of soft-parton regeneration, where enhanced q distribution significantly increases 
the p production. To learn about the effect of regeneration, we have calculated H p (x, k) with 
the soft-parton regeneration turned off, and found that for x p = 0.55 and k = 0.76 the ratio 
of the corresponding H p (x, n) values with regeneration to that without is about 2000. In 
other words without regeneration R p / P would be at the level of 2.5 x 10 -5 , which is hardly 
measurable. 

The increase of p production due to regeneration is much more than the corresponding 
increase of tt (shown in Fig. 4) for a good reason. It is not just a matter of p consisting 
of three q, while 7r having only one q. The pion recombination function is broad in the 
momentum fractions x q and x q , so with x q high it is possible for Xg to be low to reach the 
region with higher density of q. The proton recombination function is much narrower, since 
the proton mass is nearly at the threshold of the three constituent quark masses. The Xq 
values are roughly 1/3 of x p , so none of the antiquarks can have very low Xq for x p ~ 0.55, say. 
The effect of soft-parton regeneration can therefore drastically increase the p production. 

It is of interest to point out that the observed p/p ratio changes significantly with energy. 
At yfs = 200 GeV and 77 = 3.2, R p / P has been found to be 0.22, which is four times larger 
than at y/s = 62.4 GeV p2]. It implies that the degradation effect depends sensitively on 
y/s. It also means that what other ratios have been measured at y/s = 200 GeV cannot be 
used reliably as a guide for our present study at 62.4 GeV. 
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We are now able to relate the average charged multiplicity (h + + h~)/2 in the data to 
[p + p + 1.15(7r + + 7r~)]/2 that we can calculate. Since the data on the p? distribution are 
taken within the narrow band bounded by rj = 3.2 ± 0.2, we can determine the x value in 
the range of pt of interest by first identifying rj with y and use 



If we take (pr) = 1-0 GeV/c, the corresponding values of x for pion is x n = 0.4 and for 
proton, x p = 0.54, which are well inside the FR. 

The slope of the pt distribution in the semilog plot is essentially determined by the 
value of T, as prescribed by Eqs. (I2TI) and (|22|) . We find it to be T = 196 MeV. In our 
treatment here and before, the value of T incorporates the effect of radial flow and is therefore 
larger than the value appropriate for local thermal temperature that is considered in other 
approaches to recombination pm [TT]. For the values of k that can reproduce the p/p ratio 
we can calculate [p + p+ 1.15(7r + + 7r~)]/2, adjusting (pt) in fine-tuning, and obtain the two 
lines in Fig. 6. The solid line is for n = 0.76 and (pt) = 1-09 GeV/c, while the dashed line 
is for k = 0.72 and (pt) = 1.07 GeV/c. They both fit the data [3] very well. The spectrum 
being dominated by proton does not depend on k sensitively; its normalization does depend 
on the x values, which in turn depend on (pt) for fixed rapidity. For (pt) ~ 1.08 GeV/c 
the corresponding x p is ~ 0.55, which is the value we used to calculate p/p. Since the 
contributions from resonance decays have not been considered, our results for p T < 1 GeV/c 
are not reliable, and should not be taken seriously. 

In Fig. 7 we show the p/n ratio for the two cases considered above. Again, there is 
sensitive dependence on k, although not as much as in p/p. As k decreases, more soft partons 
are generated. The increase of q enhances tt production, and thus suppresses the p/n ratio. 
The dominance of proton production makes the charged hadron spectrum insensitive to the 
change in the pion sector. But the ratio manifests the pion yield directly. Currently, the 
data on the p/n ratio is still unavailable for ^ = 62.4 GeV. Since k depends sensitively on 
v^i, the ratio may be quite different from that determined at 200 GeV [Hj. To have the 
ratio exceeding 1 is a definitive signature of recombination at work. The verification of our 
results will give support to our approach of accounting for hadrons produced up to pr = 2.5 
GeV/c at rj ~ 3.2 in the absence of hard scattering. 

We note that the two lines in Fig. 7 are nearly straight, since both p and 7r distributions 
are mainly exponential, as shown in Eqs. (J2T1) and (|22l) . except for the prefactor involving pt 
for the proton. We therefore expect the measured ratio to be essentially linearly rising in Fig. 
7. But more importantly in the first place is whether the ratio exceeds 1 for pt > 1 GeV/c. 
Our concern should first be whether the regeneration of soft partons and the suppression of 
hard partons are the major aspects of physics that we have captured in this treatment. The 
precise px dependence of R p / n , i.e., whether it is linear or not, is of secondary importance 
at this point. Similarly, we expect the data to show constancy of Rp/ P for the range of px 
studied. 



x = —=e 




(27) 
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4 Conclusion 



We have extended the study of particle production in the FR and TFR to include the regen- 
eration of soft partons due to momentum degradation and to consider also the determination 
of the pt distributions. We have shown that the data of BRAHMS for forward production 
can be reproduced for all pr, when suitable values for the degradation parameter are ob- 
tained by fitting the p/p ratio. The consequence is that large p/n ratio must follow. The 
hadronization process is recombination and the pt dependence is exponential, reflecting the 
thermal origin of the partons. We predict that the exponential behavior will continue to 
higher p? even beyond the boundary separating FR and TFR. The production of protons is 
far more efficient than the production of pions. That is not surprising since it is consistent 
with the result already obtained in p] due to the scarcity of antiquarks in the FR and TFR. 
Here, the pr dependence of the proton-to-pion ratio, Rp/ n , is shown to be linearly rising 
above Pt = 1 GeV/c and can become greater than 2 above pr = 2.5 GeV/c. Any model 
based on fragmentation, whether the transverse momentum is acquired through initial-state 
interaction or hard scattering, would necessarily lead to the ratio R p / n -C 1, by virtue of 
the nature of the fragmentation functions. In contrast, Rp/ P would be around 1 if gluon 
fragmentation dominates, and be 1 if the fragmentation of valence quarks dominates, so 
Rp/ P is not the best discriminator between recombination and fragmentation. 

No shower partons are involved in the recombination process because hard partons are 
suppressed in the forward region. That is supported by the absence of power-law behavior 
in the pt dependence of the data. Without hard partons there are no jets, yet there are 
high-pr particles, which are produced by the recombination of thermal partons only. Thus 
there can be no jet structure associated with any hadron at any px- That is, for a particle 
(most likely proton) detected at, say, px = 2.5 GeV/c, and treated as a trigger particle, there 
should be no associated particles distinguishable from the background. This is a prediction 
that does not depend on particle identification, and can be checked by appropriate analysis 
of the data at hand. 
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Figure 1: The distributions of u and u in momentum fraction x with and without soft-parton 
regeneration for k = 0.6 at b = 1 fm, where k is the survival factor in momentum degradation. 
Solid (dashed-dotted) lines are for u (u) with regeneration, while dashed (dotted) lines are 
for u (u) without regeneration. 
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Figure 2: Proton distributions for (a) b = 1 fm, and (b) b = 8 fm. The solid and dashed lines 
are for the cases with soft-parton regeneration corresponding to k — 0.8 and 0.6, respectively. 
The dashed-dotted lines are for k = 0.6 without regeneration. 
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Figure 3: Proton distributions normalized by N paxt /2 for (a) k = 0.8, and (b) k = 0.6. The 
solid (dashed) lines are for the case with regeneration at b = 1 fm (8 fm). The dashed-dotted 
(dotted) lines are for the case without regeneration at b = 1 fm (8 fm) and k = 0.6. 
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Figure 4: Pion distributions for (a) b = 1 fm, and (b) b = 8 fm. The lines are as in Fig. 2. 
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Figure 5: Pion distributions normalized by N paxt /2 for (a) k = 0.8, and (b) k = 0.6. The 
lines are as in Fig. 3. 
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Figure 6: Transverse momentum distribution of charged hadrons produced in Au+Au colli- 
sions at y/s = 62 A GeV, rj = 3.2 and 0-10% centrality. Data are from [3]. Solid line is for 
case (a) k = 0.76, and dashed line for case (b) k = 0.72, calculated in the recombination 
model for the average charged hadron being approximated by [p + p + 1.15(7r + + 7r~)]/2. 
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Figure 7: Proton-to-pion ratio calculated for n = 0.76 (solid line) and k = 0.72 (dashed line) 
at i] = 3.2 in central Au+Au collisions at y/s = 62.4 GeV. 
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